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Abstract. 

We report on the controlled transport of drops of magnetic liquid, which are 
swimming on top of a non- magnetic liquid layer. A magnetic field which is rotating 
in a vertical plane creates a torque on the drop. Due to surface stresses within 
the immiscible liquid beneath, the drop is propelled forward. We measure the drop 
speed for different field amplitudes, field frequencies and drop volumes. Simplifying 
theoretical models describe the drop either as a solid sphere with a Navier slip boundary 
condition, or as a liquid half-sphere. An analytical expression for the drop speed is 
obtained which is free of any fitting parameters and is well in accordance with the 
experimental measurements. Possible microfluidic applications of the rolling drop are 
also discussed. 

PACS numbers: 47.20.Hw, 47.55.Dz, 75.50.Mm 

Submitted to: New J. Phys. 
1. Introduction 

A tiny drop of magnetic fluid responds to magnetic fields in many ways - it is a "world 
in a nutshell". Typically 1 pi of magnetic fluid (MF) contains more than 10^^ magnetic 
mono-domain particles, each with a diameter of around 10 nm, which are suspended in a 
carrier fluid like water or kerosene [1]. In the absence of an external magnetic fleld there 
is no long-range order in the MF, but when exposed to a static field the magnetic grains 
orient in part which results in a net magnetization. Application of a rotating magnetic 
field induces a torque on the suspended magnetic grains. Due to the viscous coupling of 
the particles to its surrounding carrier liquid angular momentum is transferred to the 
whole drop and an abundance of phenomena is observed. 

In a series of experiments pioneered by Bacri et al [2] a magnetic drop was levitated 
in a surrounding liquid and exposed to a field rotating in the horizontal plane. For very 
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Figure 1. Sketch of the experimental setup. For details see text. 

small angular frequencies of the field an elongated drop follows the field rotation quasi- 
adiabatically with small phase lag [3, 4, 5, 6, 7]. In the limit of high angular frequency one 
observes for small magnetic fields an oblate spheroid, for intermediate values transient 
shapes, and for large fields an oblate spheroid with "spiny starfish" appearance [2, 8, 9], 
for a review see Ref. [10]. 

Our setup, investigated in experiment and theory in this article, differs from the 
above configuration in two points fundamentally: (i) the field is rotating in a plane 
oriented vertically, (ii) the drop of ferrofluid is swimming on top of a layer of non- 
magnetic fluid. The field configuration is borrowed from a recent experiment ("the 
magnetic pump" ) where the magnetic torque drives a continuous flow of ferrofluid in an 
open duct [11, 12]. By replacing the ferrofluidic layer with a floating drop we are able 
to propel the drop with a constant translation velocity -Udrop with respect to the liquid 
surface. Moreover we could in principle manoeuvre the drop to arbitrary positions on 
the whole two dimensional liquid layer by utilizing an additional alternating field in 
y-direction. This is a new and promising technique for microfluidic applications. 

Our theoretical model describes the ferrofluid drop first as a solid sphere with a 
Navier slip boundary condition at its surface, then as a liquid (half-)sphere with own 
inner flow field. The problem is treated within Stokes approximation and the assumption 
of certain symmetries. In both cases an analytical expression for the drop speed i^dmp 
in terms of the experimentally accessible parameters is obtained. While the solution 
of the Navier slip model contains an unknown parameter, the slip length, the result of 
the liquid half-sphere model is completely free of fitting parameters and is shown to 
represent the experimentally measured dependencies very well. 

The article is organized as follows. Next we present the experimental arrangement 
together with some qualitative observations. This is followed by a comprehensive 
theoretical analysis (section 3). Subsequently the results obtained by experiment and 
theory are compared in section 4 and discussed in section 5. 
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Figure 2. The magnetic sus- 
ceptibility of the cobalt based 
magnetic liquid versus the 
external alternating magnetic 
field. The data points for the 
real and imaginary parts of 
the susceptibility arc marked 
by squares and circles, respec- 
tively. 



2. Experiment 

Our experimental setup is shown in figure 1. We place a cylindrical glass beaker in 
between a Helmholtz pair of coils that produce an externally applied field Gx{t) which 
is oriented horizontally. In addition another coil is wrapped directly around the beaker 
providing a field Gz{t) in vertical direction. Here we denote the external magnetic 
far field by G and the local one by H. A sinusoidal driving current is supplied by 
connecting the output of a function generator (Fluke PM 5138A) to one channel of a 
power amphfier (Rotel RB-1090). The input of the second channel is supplied with a 
delayed signal of the function generator. In order to allow an independent adjustment 
of both currents, a variable resistor is inserted in one driving circuit. An oscilloscope 
serves to control the phase difference of both currents. When the phase difference is set 
to 90" the coils produce a rotating field G{t) inside the beaker. Any motion of the drop 
of magnetic liquid is observed from above by means of a video camera (not shown here). 

For a good performance of the driving by the rotating field a large imaginary part of 
the susceptibility of the MF is important. Thus we have selected a magnetic fiuid based 
on air stable cobalt particles [13], which are stabihzed by oleic acid in kerosene. Figure 2 
reproduces the frequency dependence of the complex susceptibility of this fiuid measured 
by an ac-susceptometer [12]. The MF has a volume fraction of 5 % and constitutes the 
interior (i) of the drop. Its viscosity was determined to be rj^^'' = 5.4mPas by means of 
a low shear rheometer (Contraves LS40), and the density of the MF has been found to 
be p^^^ = 1.07gcm~^. 

The drop of MF has to float on top of a liquid layer of a non-magnetic fluid. The 
quantities of this fiuid outside of the drop will be marked by (o). This fiuid must not 
mix with any of the components of the MF. Moreover it must be denser than the MF. A 
per-fiuorinated hydrocarbon fiuid (Galden SV-90) proved as suitable substrate because 



Rolling ferrofluid drop on the surface of a liquid 4 




Figure 3. Drops of magnetic fluid with a volume of (a) 5 yl (see movie 1 ) and (b) 80 jil 
(mgvie2) are rolling on top of a per-fluorinated Newtonian liquid. 

of its higher density = 1.69 gcm~^, its long-term stabihty, and its non-miscibihty 
with the MF. According to the manufacturer the viscosity amounts to 77^°^ = 1.27 mPas, 
and the surface tension to 7 = 16mN/m. This fluid is poured into a cyhndrical glass 
beaker up to a height of 2 cm in order to minimize fringe effects from the bottom of the 
glass. 

At the beginning of an experiment a definite volume V of MF is put on the surface of 
the per-fiuorinated liquid with a pipette. According to the density ratio of the two liquids 
the forming drop immerses with approximately two thirds of its volume (corresponding 
to a measured maximum penetration depth of about 60 % of its diameter) . The rotating 
field generated by the coils leads to a motion of the droplets in the direction the field is 
rolling. Hence the direction of the motion can be reversed by changing the sign of the 
phase difference between the ac-fields. Under the given experimental conditions we can 
achieve droplet velocities up to a few cm/s. The good contrast between the black MF 
and the transparent hydrocarbon liquid allows an easy observation by a digital video 
camera. Two exemplary movies can be activated at figure 3. The velocity of the droplets 
was determined by extracting the time a drop takes to travel the distance of 1 cm in the 
center of the beaker. Within this distance the magnetic field varies less than 1 %. 

3. Theory 

The theoretical description of the "real" setup poses a very complicated boundary value 
problem which would have to be solved by numerical methods. In order to extract 
the essence of the effect we make some simplifying assumptions which even lead to an 
analytical solution. 

The droplet is considered to be a spherical object half-way immersed into a liquid 
with an otherwise perfectly flat surface. Effects of gravity are neglected as is the 
inertia term in the Navier-Stokes equation which is hence rendered linear. This Stokes 
approximation is in order when the Reynolds number Re is sufficiently small. Here it 
is given by Re = f2 R"^ g'^°y r]^°\ with f2 the angular velocity of the sphere, and ranges 
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Figure 4. A spherical 
ferrofluid drop with radius R 
hosts in its inner (i) a fluid 
with density g^^^ . It is covered 
from above (a) by a gas 
with density g^^\ The lower 
part of the drop is half-way 
immersed into an outer (o) 
Newtonian fluid with density 
g^°^ and dynamic viscosity 
7]^°^ The drop rotates with 
constant angular velocity O. 
The center of the sphere is the 
origin of the reference frame 
as indicated in the picture. 

between one and ten. The problem is treated within the reference frame where the 
sphere is rotating with its center at rest (cf. figure 4). In order to ensure stationarity in 
this frame, the overall forces and torques acting on the sphere must cancel out. After 
the velocity field of the surrounding liquid has been determined, its asymptotic value at 
r — > oo will give the negative translation velocity of the sphere in the laboratory frame. 

The simplest approach is treating the droplet as a solid sphere and employing 
the common no-slip boundary condition at its surface, but this would lead to a 
logarithmically divergent viscous torque [14]. It has long been shown [15], that 
hydrodynamic problems containing a moving contact line in combination with the no- 
slip condition give rise to diverging quantities due to an inherent contradiction: on the 
one hand the fluid is supposed to stick to the solid surface, and yet the line where solid, 
liquid, and gas meet shall advance on that very same surface. 

Several means have been proposed to relieve the singularity, e.g., taking into account 
a strong curvature of the fluid surface near the solid, or describing the contact region 
in terms of molecular interactions, as has been done in Ref. [16]. A straightforward 
approach is to allow a certain amount of slippage over the solid surface. As early as 
1823, a linear relation between the tangential stresses at the solid surface and the velocity 
of the latter was proposed by C.-L. Navier [17]. Although other forms of slip condition 
can be successful [18, 19] this "Navier slip" has become the most popular one and has 
since been examined and applied oftentimes. Earlier works distinguish between several 
regions where different expansions are made, and only employ the slip condition in the 
contact region itself, finally matching the solutions together [20, 21, 22]. Our treatment, 
however, will follow the lines of Ref. [23] who applied the Navier slip condition on the 
whole solid surface without seperating different regions. This is justified a posteriori by 
the fact that the slippage shows most of its impact in the direct vicinity of the contact 
line where the stresses are largest and leaves the flow field undisturbed further away, as 
will be made clear by the results of the present paper. 
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Although Ref. [23] considered a problem quite analogous to ours, i.e., the rotation 
and translation of a solid spherical object which is half-way immersed in a liquid, we 
will present the treatment in a more lucid albeit less general way that will lead to a 
closed expression for the resulting flow field which is lacking in Ref. [23]. 

The disadvantage of the Navier slip condition is that it contains a characteristic 
length Lg which is supposed to be small compared to the length scales characterizing 
the problem (in our case the sphere radius R) and essentially indicates how much the 
fluid molecules slip over the solid surface. Lg — > is equivalent to no slip, while Lg — > oo 
corresponds to completely unimpeded slip or zero tangential stress. This slip length 
does not necessarily "represent true slippage but merely recognizes the fact that the 
liquid consists of molecules of finite size" , as stated by Huh and Mason in Ref. [20] . Or 
as Cox puts it in Ref. [22]: "Slip between liquid and solid is a convenient assumption 
to get rid of the non-integrable stress singularity." Although the slip length between 
certain materials can be measured by now (see e.g. [24, 25, 26, 27]), this is of no use to 
the present problem, as the experiment does not involve a solid sphere. 

By consequence, the result of these calculations will not be entirely satisfying, so 
that a second approach is taken in which the ferrofluid drop is treated as a liquid half- 
sphere with its own inner flow fleld. In this case, the velocity flelds and also the sums 
of viscous and magnetic stresses must be continuous at the interface between the two 
liquids. Though the liquid drop cannot be described as a whole sphere but only as a half- 
sphcrc, the resultant drop speed, which no longer depends on any unknown parameters, 
represents the experimental data extremely well. This may indicate that the true flow 
field in the drop is mainly restricted to its lower part. 



3.1. Solid sphere 

The basic hydrodynamic equations are the continuity equation for incompressible fluids 
V-v = 0, (1) 
and the stationary Stokes equation 

= - Vp + TjV^v (2) 
which by eliminating the pressure can be written as 

(V X t;) = . (3) 

The velocity fleld of the non-magnetic liquid bearing the sphere is expanded in 
vector spherical harmonics according to Ref. [28, 29]. Appendix A gives the details of 
this expansion and shows how the various coefficients occurring in it are determined 
from the boundary conditions. 

When only one boundary condition is left, namely the requirement that the 
dissipating viscous torque compensate for the accelerating magnetic torque, the velocity 
components of the flow fleld below the sphere still depend on the yet unknown angular 



Rolling ferrofluid drop on the surface of a liquid 7 

velocity i? with which the sphere is rotating. The resulting expressions are (cf. Appendix 
A): 
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The flow field determines the pressure via Stokes' equation (2). Straightforward 
calculation yields 



SO that the pressure field is given by 
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3.2. Viscous torque 

The viscous torque acting on the lower half-sphere is gained from the tangential viscous 
forces 



d-Ftang = [c^r^e^ + cT^^e,^] P^d(/7 d'j? siu 1? 



(8) 
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with the tangential components of the viscous stress tensor cr^^ and a^p as defined 
in [30]. Integration over the lower half-sphere < < 7r/2, < < 27r yields the 
dimensionless viscous torque in |/-direction 
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When the doublefactorials in (10) are transformed to single factorials and Stirling's 
approximation 



£!RiA/2^/e-^ £>1 



is employed, it can be shown that the terms for large £ in the infinite series give in 
leading order 
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While is a convergent series, J^^^ •^"''^ diverges logarithmically, so here 

the necessity of the slip condition becomes manifest. 

Looking at the solution (4), (5) for the velocity field, it becomes clear that the field 
is only changed significantly near the contact line or, more generally, near the sphere 
surface: since -C R, the terms with small ^ hardly deviate from those for no-slip. 
Only when iLg/ R exceeds the order unity, the factors containing become important. 
Each term is made smaller, and the more so the greater £ becomes and, of course, the 
greater the slip length. On the other hand, the terms with large £, i.e., those which are 
influenced by the slip condition, are negligible when r ^ R. So the results with and 
without slippage would not be distinguishable far enough from the contact line. 

Figures 5 and 6 illustrate the influence of slippage in the relevant region near 
■(9 = 7r/2 for expansion orders 99 and 100, respectively. Where there is a steep descend 
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R) on 'd and therefore a large corresponding tangential 



R) for Ls 
the fluid is allowed to shp. 



0, the curves are considerably smoothed out when 



3.3. Magnetic torque 



In order to obtain an expression for the angular velocity i?, we utilize the fact that the 
viscous torque (10) must compensate for the magnetic torque which is calculated now. 
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Figure 5. Influence of 
slipping on v^'' {R) over 
i9 for = 100. Both 
the oscillations and the 
steep descent to zero are 
considerably smoothed 
out when a finite slip 
length is taken into 
account. 
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Figure 6. Influence 
of slipping on the rel- 
evant stress component 
al^\R) over for TV = 
99. The greater the 
slip length, the more are 
the oscillations damped, 
i.e., the more is the 
stress relieved. 
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The vector of the apphed magnetic field rotates within the xz-plane, generating a 
torque in y-direction, so that the external magnetic field is denoted by 

G = ^t{G}, G = G e^"*(-ie^ + e,) (14) 

with u! = 27r/ being the rotation frequency of the field and 

X = X'- iX" = Xif) (15) 

the frequency dependent magnetic susceptibility of the sphere. Concerning the 
amplitude of the magnetic field, the susceptibility is assumed to be a constant. 

The sphere is supposed to be magnetized homogeneously, having the overall 
magnetization (see for example §§8 and 29 in [31]) 



M = mc{M}, M=^e'-*(-ie, + e,) (16) 



3 



SO that the magnetic torque acting on it in the stationary state is given by [1] 

T^r., = /io^e, (M,G. - M,G.) = _ ^^^^ 

This must compensate for the viscous torque 

T^,, = -nrinR^j:{L,)ey. (18) 
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Here, the right-hand side of (10) is abbreviated by S(Ls), reminding us that it 
includes an infinite series which depends on the shp length and cannot be computed 
analytically in closed form. The equality of viscous and magnetic torques poses the 
last boundary condition which makes sure that the rotational and, consequently, also 
the translational motion of the sphere be not accelerated, and finally gives the rotation 
frequency of the sphere: 



3[(l + |)2 + (2^)2]^E(L.) -3r/S(L.) ^ ^ 

The speed with which the sphere advances on the fluid surface is given by the 
negative of the velocity field at r — >■ oo. In this limit, only the {£ = l)-tcrms remain 
and the corresponding factor from the Navier slip condition can be neglected because 
of Ls < R: 



'^drop 
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3.4- Fluid (half-)sphere 

Although a definite result has been obtained for the speed of the magnetic sphere, it 
cannot be compared to experimental data so easily. It still depends on an unknown 
parameter, the slip length Lg, which cannot simply be treated as a fit parameter. Due 
to the very weak dependence of the viscous torque on the expansion order, it poses a 
formidable numerical problem to obtain the slip length for a given torque, so it would 
be of advantage to obtain an expression for the drop speed that does not depend on 
such a parameter. 

In addition, one could expect a model containing a liquid drop to be more realistic 
than one with a solid sphere. For these reasons the ferrofiuid drop is now considered 
liquid, though still spherical, being also subject to the hydrodynamic equations like 
the surrounding liquid. The Navier slip condition is replaced by the condition of 
continuous velocities and stresses at the interface between the two liquids. All other 
boundary conditions remain as before, including the addition of the mirror image. As 
a consequence of the requirement of a fiat "surface" (f^ = at z = for all r > i?), it 
is not possible to obtain a spherical inner (i) velocity field: v^^ is rendered zero within 
the whole section z = when the corresponding outer (o) component v^^ is demanded 
to vanish on the whole contact circle z = 0, r = R. 

However, when the boundary conditions are posed in analogy to the previous 
section, a flow fleld is obtained which proves to be very useful. As the fleld becomes 
completely horizontal within the plane of symmetry, it is suggested that only the lower 
half-sphere is identifled with the ferrofluid drop, i.e., after solving the mirror image 




Figure 7. Flow field of the liquid half-sphere within the plane y = 0. 



set-up, the whole upper half-space is neglected, resulting in the flow fleld displayed in 
figure 7. 

The same differential equations (1), (3) and ansatz (A. 7), (A. 11) together with the 
requirement that the velocity be finite at r = yield for the radial functions of the inner 
velocity field {i > 0): 



fSir)^0 (21) 
filir) = qimr'+' + B,^r'-' (22) 

ht{r)=Pimr' (24) 

Starting point for the velocity components of the surrounding liquid are again the 
radial functions (A. 19) - (A. 24). For simplicity it is still assumed that the drop remains 
spherical, i.e., 

^;«(i?) =t;(°)(i?) =0 V^,yP, (25) 

instead of demanding that the normal stresses be continuous at r = i?. 

As mentioned above, the tangential components and v^p must be continuous. Due 
to the orthogonalities (A. 12) and (A. 13) this condition reduces to the radial functions 
g^rri and him being continuous. 

Furthermore, the tangential forces must cancel out at every point on the spherical 
interface so that the tangential stresses are pointwise continuous. The latter consist of 
viscous stresses cr^™^ = cr^,?/!^ and magnetic stresses [32] 

aff = lioH.H, - ^ H,H,5i, + ^ {M,H^ - M^H,) (26) 
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where i,j — x,y,z and the local magnetic field is given by 

H^^Kt{H}, ^ = ^e^-*(-ie, + e,), (27) 

3 

assuming a linear magnetization law 

M = xH. (28) 

The quantities M, G, x, and 971 are defined as in the previous section. For the 
condition of continuous tangential stresses, the symmetric part of the magnetic stress 
tensor (26) need not be considered since it is the same on both sides of the interface due 
to the usual boundary conditions for H. 

The antisymmetric part, on the other hand, is the crucial one which leads to the 
propagation of the drop. It shall be denoted by Cj-™^- Because of antisymmetry in 
addition to Hy — My — 0, only one independent cartesian component is left: 

= ^ , = -Tl (29) 

This gives the tangential magnetic stresses 

alf = mcos(p (30) 
(jM = -9Jtcos^?sin(^. (31) 

Now the boundary condition reads 

Fi7l{R)=Fip{R) + Fj;';;-\R) (32) 

because the accelerating magnetic force must be compensated by the viscous ones. 
With F^/^ = Uy^i^ n-e,. and the convention that the surface normal n = +6^ for a force 
that acts on the outer surface and n = — for a force that acts on the inner surface, 
this yields in terms of stresses 

^i^. + '-i?;;^(^)-^i^n^) = (33) 

for all 1?, Lp. As before, the viscous force in x-direction must vanish. The resulting 
expressions of the components of inner and outer fiow field are given explicitely in 
Appendix B. 

Again, the speed of the drop in the laboratory frame is obtained by evaluating the 
negative of the outer velocity field at r — > oo, giving 

^SL--- n e.. (34) 

drop 2 2r/(°) + 3r/(') ^ ' 

Although this result looks very similar to the one obtained in the previous section, 

, _ 4 TIR 

it clearly has two advantages. First, it purely consists of parameters that are 
experimentally measurable or tunable (sphere radius R, viscosities rj, and via 9Jl 
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susceptibility x external magnetic field amplitude G). Second, there is no need 
of calculating numerically an infinite sum. 

Since no singularity has occurred in the scope of the calculations for the liquid 
sphere, it can be compared to a model where slipping is taken into account. The stresses 
which diverge within the framework of the very rigid no-slip condition are relieved both 
when the surrounding fluid is allowed to slip over the solid and when the solid is replaced 
by an elastic or, as in our case, viscous medium. Indeed, the crucial viscous stress 
component o"'"^ = cr^^^'°^ is essentially identical to the one obtained from the velocity 
field with Navier slip, the only differences being constant factors, at least when £^ 1: 
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4. Comparison of experimental and theoretical results 

The main result of the theory for the drop speed (34) reads exphcitly 

liq ^ _R I^qG'^ X^^ (oo\ 

^drop - 4 2rj(o) + 3^(i) • (1 + 1)2 + (2|^)2 • ^-^^^ 

It can be well compared with the data obtained in the experiments. They have 
been measured following the procedure described in section 2. Figure 8 presents a plot 
of the drop velocity versus the magnetic field amplitude G for a driving frequency of 
/ = 0.8kHz. We have put a droplet of volume V — corresponding to a sphere 
of radius R ~ 1.1mm, on top of the liquid layer. The measured velocities (marked 
by full circles) show a monotonous increase with G. The solid line gives the values 
of (38), taking into account the viscosities of the ferrofluid, i]^^^ = 5.4mPas and of 
the liquid below, which amounts to 1]^°^ = 1.27 mPas. The driving frequency enters 
into expression (38) only via the real and imaginary part of the magnetic susceptibility 
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Figure 8. The drop 
speed in dependence of 
the magnetic field ampli- 
tude G for / = 0.8 kHz 
and V = Syl. The 
blue circles mark the 
measured data, the red 
line gives the theoretical 
curve according to (38), 
taking into account the 
proper material values. 




Figure 9. Drop velocity 
versus drop radius for 
an alternating magnetic 
far field with G = 
0.844 kA/m and / = 
0.8 kHz. The blue dots 
mark the experimental 
results, the solid line the 
theoretical outcome. 




0.5 1 1.5 
drop radius R / mm 



which were determined as = 4.66 and x'' = 3.25, respectively, for the given frequency 
(cf. figure 2). As can be seen, the values for the liquid half-drop solution represent the 
given experimental data extremely well. 

In a series of measurements different drops with a volume ranging from 1 to 50 pi 
were investigated. For comparison with theory we assume a spherical symmetry and 
estimate the drop radius R from the drop volume V. As shown in figure 9, the measured 
drop velocity increases with the radius of the drops. The solid line marks the result of 
(38) for an amplitude of G = 0.884 kA/m, as set in the experiment. Again we find a 
quantitative agreement of the half-drop solution with the experimental data. 

As a further parameter the driving frequency / was varied in the experiment. When 
the frequency dependence of the drop speed was determined, the vertical field was fixed 
at Gz = 0.844 kA/m. However, the frequency dependent inductance of the outer coils 
did not permit to keep at this value for the whole range of frequencies (the ratio 
Gx/Gz is indicated at the r.h.s. of figure 10). In order to obtain a magnitude which is 
independent of G, we introduce the reduced velocity 

"^'^''^Pi^/zoG.G,' ^^^^ 
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X 10^ 




frequency f / kHz 



Figure 10. Frequency dependence of the reduced drop velocity u for V ~ 5 pi. The 
full circles mark the experimental data, the solid line gives the theoretical curve where 
the measured frequency dependence of x"{f) has been plugged in. For all 

data Gz was fixed to 0.844 kA/m, but Gx was decreasing with increasing /. The green 
open squares are indicating the actual ratio Gx/Gz. 

where Gx denotes the horizontal and Gz the vertical field amplitude. Within the 
linear regime this quantity should be independent of the choice of the amplitudes. 
Figure 10 shows an increase of the reduced drop velocity (marked by solid circles) 
up to a maximum at / = 10 kHz. The theoretical values (solid line) stem from (38), 
where the material parameters and the measured frequency dependence of the complex 
susceptibility x'(/) +^X''(/)) as presented in figure 2, have been utilized. In order to be 
able to compare the predictions with the experimental results, v^^^^ is scaled according 
to (39). We observe a good agreement up to a frequency of about / = 1.5 kHz. Beyond 
that point, the theoretical curve deviates from the experimental results. The former 
shows a maximum at about / = 3.5 kHz, while the measured velocity is largest at 
/ = 10 kHz, and the maximum values differ by a factor of two. 

5. Discussion and Conclusion 

The measured propagation velocity of the droplet shows a parabohc dependence from the 
magnetic field amplitude, and a linear dependence from the radius of the droplet. Both 
experimental observations are quantitatively described by the liquid half-drop solution, 
without any free fitting parameter. The theory just needs the magnetic field amplitude, 
the complex susceptibility and the viscosities of both fiuids {i.e., the ferrofiuid and the 
liquid layer). Taking into account the over-simplifying assumption of a half-spherical 
drop, the theory describes the experimental data remarkably well for driving frequencies 
up to 1.5 kHz. 

For higher driving frequencies, however, (cf. figure 10) a discrepancy between 
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experiment and theory of up to 100 % is observed. This discrepancy may have several 
origins. Firstly, due to experimental characteristics, the rotating magnetic field becomes 
elliptical. Following [9], the nonlinear effects of an elliptical field are expected to diminish 
the flow within the droplet. This, however, does not explain our experimental data, 
which overcome the predictions by theory. Of course our experimental situation differs 
from that of [9] where an elliptical drop can freely rotate in the horizontal plane. In 
our case the horizontal surface is pinning a free rotation of an elliptical droplet in the 
vertical plane. 

Secondly, for higher driving frequencies the liquid-liquid interface of a fully 
immersed drop develops spikes and resembles a "spiny starfish", as reported in 
Refs. [2, 9]. This may also happen for the lower part of our half-immersed, swimming 
drop. A complex interface of the two liquids may enhance the interaction in between the 
fluids and thus increase the propulsion - similar to a paddle wheel of a Mississippi steam 
boat. This can of course not be covered by the simplifying model ansatz. The shape 
and dynamics of the liquid-liquid interface shall be studied in forthcoming experiments. 

The main achievement of the article is that rotating flelds can transport ferrofluidic 
drops. Our experimental results can be quantitatively explained without any free fltting 
parameters. 

Moreover the theory gives an explicit solution of the flow flelds both for a rotating 
solid magnetic sphere and a spherical ferrofluid drop which both arc half-way immersed 
in a liquid. The similarity of the final results of both cases demonstrates the equivalence 
of Navier slip at a solid surface on the one hand and the continuity of tangential stresses 
at a fiuid-fiuid boundary on the other hand. 

For a quantitative description of " magnetic pumping" by means of a rotating field 
a droplet is more suitable than a plain ferrofluidic layer [12]. For the droplet one does 
not need any tracer particles (the droplet is its own tracer), and the demagnetization 
factor of an eUiptical droplet is well deflned. 

Future experiments shall unveil whether the half-drop model works also in the pico- 
litcr range. Here the dimensioning of droplets is very precise (see e.g. Ref. [35]) and 
their position may be detected by magnetic sensors [36]. Taking advantage of (38) one 
may even select the size of the generated droplets by their speed. 

We propose that the controlled transport of small amounts of liquid to any desired 
position on top of a liquid two dimensional layer is a promising technique for microfluidic 
applications. There ferrofluidic drops are commonly manipulated utilizing local fleld 
gradients, which are locally created by embedded wires [33] or planar coils [34]. In 
contrast, our driving technique yields a constant drop velocity globally, i.e. on the 
complete surface. 
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Appendix A. Explicit computation of the flow field below the solid sphere 
The velocity field is expanded in vector spherical harmonics according to [28, 29] 



v{r, ^:^) = Y,^ { ^rfem{r) + gim{r)rV + him{r)r x V Wem{^, ip) (A.l) 

with the normalized spherical harmonics Y^^ and the Legendre functions as 
defined in [37] for £ > and < m < £: 



YU^,^) = X e'"^^P.„.(cos^) (A.2) 

= i^^„e^"^^P^^(cos^?) (A.3) 

y,,_U^,(^) = (-iry4(^,(^) (A.4) 

PUcos^) - ^-^i^^^r^^^^^isin^r (A.5) 

When the expansion (A.l) is put into (1) and (3), these partial differential equations 
for the vector v are transformed to ordinary differential equations for the scalar radial 

functions fim, gim, and him- Before this is done, equation (A.l) is simplified by several 
means. 

With the Nabla operator in spherical coordinates 

V = erdr + -e^d^-\ e,. (A. 6) 

r rsmi7 

where dj = d/dj, the velocity components read: 

oo +e 

^^r- = XI XI f^rnYem (A.7) 

e=o m=-e 

\ \ ^ \ ^ I / 1,1 ■ n-'-im 




+ h,ra sin^ I ) (A.8) 

With (A.4) and the fact that the velocity field is real valued it follows 

g,,-m = {-iTgem. he,-m = {-irhem ■ (A.9) 
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Furthermore, when the symmetry of the problem with respect to the xz-plane, i. e. 

v^{-f) = v^{ip), v^{-(p) = -v^{ip) (A.IO) 
is taken into account, it can be shown with the aid of relations (A. 4) and (A. 9) that 



oo i 



/ , / , '^Kim { 9em I . , , -m 



m 



v. 



sm{m(p)- — -Pem 
smv 



h 



£m 



cos{m(p)- — -Pem 
smv 



sm{m(p) d^P(,m 

= 2 ^ ^ ^gimAlm + himBirn^ . (A.H) 

l=Q m=-e 

The prime at the second sum indicates that the terms with m = are divided by 

two. 

When the boundary conditions are applied it will be important that the two velocity 
components of (A. 11) always be considered together, because A^^. — -^imi"^,^) and 
— Bim{'&, f) fulfil the orthogonality relations 



{Aim, Bi'rn') = 

{Aim, Ae'rn') = {Bem, Bi'm') = 7- 

with the vector inner product 

(Xi,X2):=/ dip d^sm^{X*fX2, 



(A.12) 
(A.13) 



(A. 14) 
Jo Jo 

where * denotes the complex conjugate and ^ the transpose of the vector. By 
computing the inner product of A^'m' or .B^/^/ with (A. 11) one can reduce the infinite 
series to one function ge'm'{T^) or he/m'ir), respectively. If and were considered 
seperately, it would not be possible to get at the radial functions, because d^P^'m' a-nd 
i^T^P^OT alone are not orthogonal. 

Now putting the expansions (A. 7) and (A. 11) into the basic equations (1) and (3) 
gives the following ordinary differential equations for the radial functions with £ > 
{goo{R) = hoo{R) = can be assumed w.l.o.g.): 



/oo + ~/oo — 
r 



' f"" I " f" _L _ 

= 



(A.15) 



6 



- 1 



+ 1) 



£{£ + !) 

''"fern + 
+ 1 



fim o ft 



2 J im 



+ 1) -2 



fi 



em 



tm 



-h 



£m 







(A.16) 
(A.17) 

(A.18) 
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These equations are solved by a power law ansatz which together with the 
requirement that the velocity be finite as r — > oo leads to 

himir) = e>0 (A.19) 

/oo(r) = ^ (A.20) 

fimir) ^b,m + — + ^ (A.21) 





1 '^Im 


r 






dim 




2r3 



^i.(r) = + ^ - ^ (A.22) 
and for i > 1: 



fem{r) 
gtm{r) 



+ 



(A.23) 
(A.24) 



£(£+1) 

The coefficients a^^, q^, and are determined by successively applying the 
remaining boundary conditions. In the following section, the ferrofluid drop will be 
treated as a solid sphere. Its angular velocity J? is introduced as a parameter that will 
have to be determined by the equality of viscous and magnetic torques generated by 
external field and surrounding liquid. 

It should be noted here, that the orthogonality relations (A. 12) and (A. 13) would 
not be valid if the i^-integral within the scalar product (A. 14) were only carried out 
up to ■»? = 7r/2. On the other hand, the hquid only occupies the lower half-space in 
the given problem, so we perform a little trick in order to be able to integrate over 
the whole sphere, i.e., we take advantage of our equations being linear and employ the 
superposition principle by adding the mirror image of our problem with respect to the 
xy-plane (fluid above, void below the sphere). The problem can be solved in this way 
and the resulting flow field in the upper half space is simply neglected in the end. 

Within the framework of this "mirror image construction" the following boundary 
conditions are employed: 

• Navier slip at the sphere surface 



r=R 



1 



vAr 



.v^[r 



U 



(A.25) 



with the shp length Lg <^ R and the velocity U of the sphere surface 




implying 



U 



v^{r = i?) = 0. 



for 1? = 7r/2 
for ^ < 7r/2 
for > 7r/2 



(A.26) 



(A.27) 
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(A.28) 



/'27r P'^/^ 

Fi = d(p d'd sin i? arj (r = R) ej ■ = 
Jo Jo 



(A.29) 



with i e {x,y,z}, j G {r,{},ip}, and e^, the unit vectors in respective direction. 
The relevant components of the viscous stress tensor are taken as defined in 
[30]. 

As is obvious from the given symmetry, only F^. will be different from zero and 
thereby determine the last coefficient. 

Since the magnetic field only creates a torque but no linear force, this boundary 
condition provides the requirement of unaccelerated translational motion. 



Appendix A.l. Applying the boundary conditions 
bcapp 

The first coefficients are determined by the r-component of the Navier slip condition 



oo +1 



=0 m= 



and the orthogonality of the scalar spherical harmonics Y^m 



(A.30) 



fUR)= V£,m 
doo = 

^ dim — ^R^'blm — R^Cim 
dim — —R^Cim, i > 1 



(A.31) 
(A.32) 



The coefficients and are obtained by applying the appropriate vector inner 
product to the and (/^-component of the Navier slip condition 



1 + ^- L^dr 



e=i m=0 



em 



(cos{m(p) 
—m 
sm mv?^— 
sm'jy 





/?i?(cos</?e^ — cos sin (/? e^) 
— i?/2(cos</7e^ — cos sin (/? e,^) 



sin(m(/?) 
for ^ = n/2 
for ^9 < 71 /2 
for 1? > tt/2 



Im 



R 



(A.33) 
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which is done here exemplary for the scalar product with A^r^, (i?, (p) as defined in 
A. 14. The orthogonalities of the sine and cosine functions yield 



1 + § ) 9imiR) - U gUR) = V m 7^ dzi 



and 



7r/2 

'o 



J-k/2 

Now from [38] one finds 

d^Pa + cot mi = ^(^ + l)P^o 

and 

^ duPt{u) = j duPi{u), 



+ 1) = 
a^ + coti^jPalcosi?) 
a^ + coti?lPa(cosi?). 



^ even 
£odd 

so that with the definition of K^^n according to (A. 3) one obtains 

Ml even 



^ duP^{u) = - j dup,{u) = (-i)^|jq:Yy!!' 



l + ^)5a(it:)-i^s^a(it:) 



and for odd £ 



(2/ + 1)^ , , m (£-2)!! 



+ 1) 

With (A.22) and (A.24) this gives in detail 



+ 1)! 



^10 



hi 



,±1 



3 M 
- + 3— 
2 R 

3 Ls 
- + 3— 
2 R 



+ 



ClO 

R 



1 + 3 







+ 



Cl,±l 



Qm = 

Q,±i = 



Vm ^ ±1 
V i even 



1 + 3 



37r 




(A.34) 



(A.35) 



(A.36) 



(A.37) 
(A.38) 

(A.39) 

(A.40) 



(A.41) 
(A.42) 

(A.43) 



l + (2£ + l)^ (£ + 1)!!' 



£odd. 
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gioir)Kio a^Pio(cosi?) _^ + XI X] 9imir)Kim d^Pimicos'd) 

~2 m=±l t=i 
I odd 

oo 

+ X X mhim{r)KimPtm{^) = 0. 



(A.44) 



m=±l 1=2 
t even 



The sums vanish completely due to properties of the Legendre functions at zero 
[38] , so that only the flrst term remains, giving 



310 



1 + 



1^ 



+ 



ClO 

2r 



1 + 



0. 



Since this equation must be valid for arbitrary r it follows 6io = = cio- 
In order to evaluate the force condition 



Ft. — 



2-K j-K/2 

dip I d'j^sin'j? 



(Trr{R) Sin'J?C0S</7 



+(Tri?(-R) COS 'd COS if — ar<p{R) sin If 

the following integrals are needed: 



d-d sini? 



1 



sin'j? 



+ cos "d 



Then the last coefficients are given by 



bi,±i = T 



QR 



1 + 2^ V 6 



Appendix B. Resulting flow fields for the liquid half-sphere model 
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